Analog of photon-assisted tunneling in a Bose— Einstein condensate 
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We study many-body tunneling of a small Bose-Einstein condensate in a periodically modulated, 
tilted double-well potential. Periodic modulation of the trapping potential leads to an analog of 
photon-assisted tunneling, with distinct signatures of the interparticle interaction visible in the 
amount of particles transferred from one well to the other. In particular, under experimentally 
accessible conditions there exist well-developed half-integer Shapiro-like resonances. 

PACS numbers: 03.75.Lm, 74.50.+r 



In a recent experiment, it has become possible to 
monitor the Josephson-like oscillation of a sample of 
about 1000 Bose-Einstein-condensed atoms in an opti- 
cal double-well potential, and to observe in situ both 
the evolution of the atomic densities, and of the relative 
phase between the condensates in both wells Q. In the 
present Letter, we propose to extend such experiments by 
modulating the trapping potential periodically in time. 
As we will demonstrate, such modulated double- well con- 
densates exhibit certain parallels to photon-assisted tun- 
neling in microwave-driven superconducting Josephson 
junctions Q, although the "photons" now tend to fall 
into the lower Kilohertz regime. 

We consider a symmetric double-well trap filled with 
N Bose particles at temperature T = 0, as realized ap- 
proximately in the Heidelberg experiment Q]. Denot- 
ing the tunneling splitting between the lowest pair of 
single-particle energy eigenstates by HfL and adopting 
the common two- mode approximation |3j, LJ, La] , the sys- 
tem's Hamiltonian takes the idealized form 

2 ' (1) 

where a™ is the annihilation (creation) operator for a 
Boson in the i-th well, satisfying the commutation rela- 
tion [a», at] = Sij = 1,2). The on-site interaction 
energy of a single pair of Bosons occupying the same well 
is 2hn, proportional to the s-wave scattering length of 
the particle species. This deceptively simple model (Q, 
describing an unforced Bosonic Josephson junction, has 
recently been studied in considerable detail 0, Q ■ The 
key parameter governing its dynamics is the dimension- 
less ratio Nk/SI. While the well-understood mean- field 
approximation 0, 0, Hcj corresponds to the limit N — > 00 
and k — > 0, taken such that the product Nk remains con- 
stant, present experiments start to explore the dynamics 
beyond the mean- field regime 0. 

We extend this system by introducing a twofold bias: 
We assume that the two wells are tilted such that their 
bottoms are misaligned in energy by an amount 2h/iQ, 
which should be sufficiently large compared to the tunnel 
splitting Ml, so that the usual Josephson oscillations are 
strongly suppressed. In addition, we propose to modu- 



late this tilted double well periodically in time, such that 
the individual wells are shifted sinusoidally up and down, 
in phase opposition to each other and with angular fre- 
quency u>, by an amount h[i\. Such a modulation of an 
optical double well can be achieved by periodically shift- 
ing the focus of a blue-detuned laser, which creates the 
barrier between the two wells, with the help of a piezo- 
actuated mirror 0. Thus, we are led to the explicitly 
time-dependent many-body Hamiltonian 

H(t) = Hq + fi(/io + Ml sincji) (a^ai — a\aij . (2) 

As follows from elementary estimates [jj 0, EI , the un- 
perturbed two-mode Hamiltonian can be trusted as 
long as the particle number N does not exceed the ratio 
of the characteristic length scale L of one of the wells and 
the scattering length a; for /mi-sized traps filled with al- 
kali atoms, this ratio L/a will typically be on the order 
of 1000. However, the ability to reduce atomic scattering 
lengths by means of Feshbach resonances 0] allows one 
to extend the two-mode approximation to much larger 
samples. Moreover, both Rabi frequencies /irj and /ii and 
the modulation frequency w should remain restricted to 
a few times the tunneling frequency 17. With tunneling 
times on the order of some 100 ms 0,0], this leads to 
modulation frequencies in the low kHz-regime. 

In the following, we investigate the dynamics of the 
model J2J) • In our calculations we assume that initially, 
at time t = 0, all N Bose particles are prepared in the 
lower well (that is, in the well tilted downward by — ft/io). 
When solving the time-dependent Schrodinger equation 
numerically, we record the expectation value 

(J z )(t) = i<^(t)|o^Oi - a\a 2 \^{t)) , (3) 

which quantifies the imbalance of the numbers of parti- 
cles found in the individual wells in the course of time, 
such that (J z )(t )/N = +1/2 (-1/2) means that all par- 
ticles occupy the lower (higher) well at i . We keep the 
static tilt fixed at three times the tunneling splitting, 
2/io = 317. As shown in Fig.^ this value allows merely a 
small-amplitude oscillation of the population imbalance 
when there is no periodic modulation (dashed line). 
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FIG. 1: Time evolution of the population imbalance {J z )/N 
for a Bosonic Josephson junction with N = 1000 particles, 
interaction parameter Nk/Q — 0.5, and tilt 2fj,o/fl = 3.0. 
The dashed line refers to the system without periodic forcing; 
the two full lines are obtained for periodic modulation with 
one-photon resonant frequency u) /Q — 3.0. The strongly oscil- 
lating line corresponds to the driving amplitude 2/xi/w = 1.8, 
close to the first maximum of the Bessel function Ji (inset: 
solution to the nonlinear Schrodinger equation for this situ- 
ation), while the full line at the top of the figure results for 
2/ii/oj = the first zero of Ji. 

In order to condense information about the system's 
response to bias with various parameters in a single view- 
graph, we characterize each quantum trajectory by its 
time-averaged imbalance 

(J,)t=±f ^(J,)®, (4) 

employing the averaging interval At — 100/57. Figure[21 
depicts results of such calculations for a periodic mod- 
ulation with fixed scaled amplitude 2\x\/lo — 0.5, and 
frequencies to ranging from 0.5 O to &Q. By construc- 
tion, values of (J z )t/N close to 0.5 indicate that the con- 
densate remains trapped almost entirely in the initially 
occupied well during At, whereas values close to zero 
signal that it visits both wells about equally. From the 
physics of superconducting Josephson junctions, one ex- 
pects photon-assisted tunneling to take place when the 
energy of an integer number of photons matches the static 
tilt, nhuj — 2hn , as in the case of Shapiro resonances 0- 
Figure |3 indeed reveals a pronounced single-photon res- 
onance at lu/Q w 3.0, clear signatures of two-photon- 
assisted tunneling at co/Q « 1.5, and even traces of a 
three-photon process at oj/Q « 1.0. However, the ac- 
tual resonance condition does not refer directly to the tilt 
2h/j,o, but rather to the tunneling frequency of the con- 
densate without periodic forcing: Within the mean-field 
approximation, the dynamical system (J2J corresponds to 
a driven nonlinear pendulum ||; in such a system, a res- 
onance can occur when a rationale multiple of the mod- 
ulation frequency uj matches the oscillation frequency u>q 
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FIG. 2: Time-averaged population imbalance {J z )t/N for the 
driven Bosonic Josephson junction with N = 1000 particles, 
tilt 2//o/fi = 3.0, and scaled driving amplitude 2/ii/w = 0.5. 
The interaction parameters are Nk/Q — 0.1 (full line), 0.3 
(short dashes), 0.5 (long dashes), and 0.7 (dots); the averaging 
interval is A(Qt) — 100. The inset highlights the half- integer 
resonance at uj/2 ~ uq(Nk/Q,). 

of the unforced pendulum [l3j . This resonance condition 
obviously is effective even when there is no atom-atom 
interaction. For Nk/Q = 0, one has uj (0) 2 — 4^q + Q 2 , 
shifting the one-photon resonance from the Shapiro value 
uj/Q = 2/j,q/Q = 3 to m/Q ~ 3.162. Even with weak in- 
teraction, Nk/Q = 0.1, this shift remains visible in the 
location of the main resonance peak in Fig. 

With increasing interaction Nk/Q, the Shapiro- like 
resonances depicted in Fig. shift downwards in fre- 
quency. (They shift upwards, by different amounts, if 
initially the higher well is occupied.) The mean-field 
approximation allows one to explain this observation 
quantitatively by expressing the oscillation frequency 
loq(Nk/Q) of the undriven, interacting system in terms 
of elliptic integrals We have checked that the con- 
dition nto m loq(Nk/Q) with n = 1, 2, and 3 reproduces 
the locations of the one-, two-, and three- photon-like res- 
onances quite well. 

Figure |3 also reveals a well-developed half-integer res- 
onance at lu/2 « u)q(Nk/Q), corresponding to n = 1/2. 
Such subharmonic resonances can again be understood 
with the help of the nonlinear Schrodinger equation, 
which describes the dynamics of the driven condensate 
within the mean-field approximation: The nonlinearity 
effectively generates both harmonics and subharmonics of 
the driving frequency. For a matter-of-principle demon- 
stration, Fig. |3 shows the time-averaged imbalance pro- 
vided by the nonlinear Schrodinger equation for stronger 
forcing with scaled amplitude 2/xi/w = 1.8 and rather 
weak interaction Nk/Q = 0.1, together with the corre- 
sponding data for a non-interacting system subjected to 
the same forcing. While the mean-field result exhibits 
half-integer resonances with n — 3/2 and n = 5/2, these 
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FIG. 3: Time-averaged population imbalance {J z )t/N for the 
driven Bosonic Josephson junction with tilt 2fj,o/fl = 3.0 and 
scaled driving amplitude 2pi/oj = 1.8. The full line is ob- 
tained from the nonlinear Schrodinger equation for Nk/Q = 
0.1; the long-dashed line refers to a non-interacting system. 
Also shown is the result for a system with only N = 2 par- 
ticles, keeping Nk/Q — 0.1 fixed (short dashes; displaced 
downwards by 0.1 for clarity). 



resonances vanish when Nk/CI = 0. On the other hand, 
the graph of the imbalance 10} provided by the Af-particle 
Schrodinger equation for N = 1000 (not shown) is prac- 
tically indistinguishable from that of the mean-field data 
displayed in Fig. [31 and even a two-particle system yields 
a pronounced peak for n = 3/2 when Nk/Q is adjusted to 
0.1. In general, while the time-dependent A^-particle im- 
balance J2J follows the mean-field data only in the limit 
N — > oo (with k adjusted such that Nk stays constant), 
but soon deviates from its mean-field counterpart when 
N is of the order of 1000 (cf. Fig. GJ, it is remarkable 
that its average (@J comes fairly close to the mean-field 
approximation even when TV is much smaller. 

A powerful theoretical tool for studying the many-body 
dynamics of the driven Bosonic Josephson junction is 
provided by quantum Floquet theory: Since the Hamil- 
tonian (0) is periodic in time, H(t) = H(t + T) with 
period T = 2h/oj, there exists a complete set of solu- 
tions to the time-dependent Schrodinger equation of the 
form \ip m (t)) — \u m (t)) exp(— ie m t/K) , with the Floquet 
functions \u m (t)) = kuj(t+T)) sharing the T-periodicity 
of the Hamiltonian |14| . The quantities e m , which serve 
to describe the time evolution in a similar manner as 
energy eigenvalues do in the case of energy eigenstates, 
are called quasienergies. The existence of these time- 
dependent many-body states, which incorporate both the 
interparticle interactions and the external forcing in a 
non-perturbative manner, is not tied to the particular 
model Hamiltonian J5J, but stems solely from the tem- 
poral periodicity of the forcing, exactly as the existence 
of Bloch states in solid state physics is solely due to the 



While the exact calculation of Floquet states and 
their quasienergies in general has to ressort to numerical 
means 11], the particular case of an integer resonance, 
nuj = 2/io, allows for an instructive analytical approxi- 
mation. Provided the static tilt is sufficiently large to 
substantially reduce the unperturbed Josephson oscilla- 
tions, and assuming low resonance order n, so that cj/Q. 
remains large compared to Nk/Q, the guasienergies of 
the driven junction (J2J) are approximately given by the 
energy eigenvalues of the undriven junction but with 
an effective tunneling frequency 



^eff — QJn { — — - 



(5) 



where J n {z) is an ordinary Bessel function of order n. 
This renormalization of the tunneling frequency under 
resonant driving resembles the renormalization of atomic 
^-factors by radiofrequency fields ^(|, or the coherent 
destruction of single-particle tunneling [l7| . 

To demonstrate that this concept is meaningful under 
experimentally accessible conditions, we haved displayed 
in Fig. H the evolution of the population imbalance © 
under the influence of a strong modulation: The lower, 
strongly oscillating full line results from one-photon res- 
onant driving with scaled amplitude 2fii/uj — 1.8, which 
is close to the first maximum of J\. Correspondingly, 
one observes photon-assisted tunneling with almost com- 
plete population exchange during the initial stage, level- 
ing down to smaller-amplitude oscillations which average 
to zero. As remarked above, the solution to the nonlinear 
Schrodinger equation for the same situation (inset) devi- 
ates strongly from the iV-particle solution after only a few 
periods. The upper full line in Fig. ^ results from one- 
photon resonant driving with 2/Xi/w ~ 3.832, the first 
zero of J\\ Since the effective tunneling frequency (jSJ) 
vanishes then, the particles remain caught almost per- 
fectly in the initially occupied well. 

The renormalization Q prompts at an interplay be- 
tween photon-assisted tunneling and an effect known as 
self-trapping. It has been established that when the ab- 
solute value of the interaction parameter \Nk/Q\ exceeds 
a certain critical value, which equals 1 for maximum ini- 
tial population imbalance, the nonlinear Josephson os- 
cillations occuring in the undriven junction become sup- 
pressed sua ; this self-trapping phenomenon has by now 
been observed in the pioneering experiment [l|. There- 
fore, employing driving frequencies sufficiently high for 
the modulation-induced renormalization © to hold, and 
starting from an initial state with all particles in one well, 
we expect self-trapping to occur in a resonantly driven 
Bosonic Josephson junction when 

Nk 



spatial periodicity of a crystalline lattice 



defy 
15]. 



fiJ n (2Mi/w) >1 - (6) 

This expectation is confirmed in Fig. 0] where we have 
plotted time-averaged population imbalances as func- 
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FIG. 4: Time-averaged population imbalance {J z )t/N for the 
driven Bosonic Josephson junction with N = 100 particles, 
tilt 2^io /Q = 3.0, and one-photon resonant driving frequency 
ui/Q = 3.0. The interaction parameters are as in Fig. [5] 



breathing JSJ) of the well-to-well coupling under strong 
resonant driving. However, driven Bosonic Josephson 
junctions offer further peculiarities of their own, such as 
interaction-induced resonance shifts, and clearly devel- 
oped half-integer Shapiro- like resonances, which should 
be verifiable under presently accessible laboratory con- 
ditions. Subharmonic Shapiro steps in microwave-driven 
superconducting point contacts have recently been as- 
cribed to multiple Andreev reflections . The fact that 
related phenomena occur in modulated optical Bosonic 
junctions, where all of the parameters N, fl, k, fXo, Hi, 
and lo can be tuned separately, opens up far-reaching new 
avenues of investigation. 
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tions of the scaled driving amplitude, again for the exam- 
ple of a one-photon resonance. The number of particles 
has been reduced to N = 100 here: As far as the aver- 
ages @ are concerned, this yields practically the same re- 
sults as N — 1000, and even as the nonlinear Schrodinger 
equation, despite the differences illustrated in Fig.^ For 
very low driving amplitudes, 2/ii/a; < 0.1, the forcing is 
hardly effective, so that the static tilt keeps the parti- 
cles in the initially occupied well. With increasing am- 
plitude, photon- assisted tunneling sets in. For 2/xi/u; 
larger than about 1.5, the four curves drawn in Fig. 01 
behave rather differently: For low interaction parame- 
ter Nk/CI = 0.1 (full line), the trapping condition © 
is effective only in comparatively small neighborhoods of 
the zeros j\ t \ ~ 3.832, jx,i — 7.016, and ji^ ~ 10.173 
of J%. When Nn/Q is increased to 0.3 (short dashes), 
these neighborhoods become wider; in the entire interval 
between and j'1,3 the condition © is met, so that 
the imbalance does not approach zero. For Nk/H = 0.5 
(long dashes) self-trapping occurs already between the 
first and second zero, while Nk/H = 0.7 (dots) does not 
admit full population exchange even for scaled ampli- 
tudes below jx,\i in agreement with the inequality 
Thus, the interactions among the Bose particles exert a 
strong influence on photon-assisted tunneling, leading to 
a distinctly non-monotonous dependence of the popula- 
tion transfer on the driving amplitude, and allowing sub- 
stantial transfer only in limited subsets of the parameter 
space. 

In conclusion, we have suggested to extend the on- 
going experimental investigation of condensate tunneling 
in optical double wells [l| by a time-periodic modula- 
tion of the trapping potential, thereby opening up paral- 
lels to the ac Josephson effect. Many-body Floquet the- 
ory leads to the prediction of the amplitude-dependent 
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